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1 Introduction 



In several recent papers [l]-[6] the hypothesis of a universal low-energy behaviour of the QCD 
running coupling has been used to estimate power-behaved corrections to various QCD observ- 
ables. While the only essential parameters are a few moments of the coupling in the infrared 
region, the absence of a satisfactory explicit model for the non-perturbative contributions has 
made the arguments hard to follow and the approximation scheme unclear. In the present paper, 
a simple model is introduced with sufficiently good properties to allow quantitative comparisons 
of various proposed schemes for the estimation of power corrections. 

Let us start by recalling that the running coupling has the general form 

as{k^) ^ ^a{k^/A^) (1.1) 



where for Nj active flavours 



and to 1-loop order at high energies 
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(1.2) 



a(x) ~ ^ . (1.3) 

A fundamental problem with simply continuing this expression to low energies is the Landau 
pole at k'^ = Is? {x = \). On general causality grounds, at least if the running coupling is to be 
related analytically to some physical observable, we expect singularities of the function as{k'^) 
to occur on the negative real axis only. 

A number of explicit models have been proposed [7]- [10] for including non-perturbative 
contributions at low k'^ and cancelling the Landau pole. A particularly simple one with many 
good features is the 'analytic model' [8, 9], 

a{x) = J- + . (1.4) 

In a; 1 — x 

This expression has no Landau pole, only a branch point at k'^ = 0. However, the extra term 
introduces a 1/A;^ correction to the logarithmic running of the coupling at large k'^. This would 
have to be cancelled somehow in quantities that are proportional to ctgiQ'^) in leading order but 
are not expected to have l/Q"^ corrections, such as the e+e~ total cross section. In addition 
the numerical value of ag becomes rather large, as(0) = l/6o = 1-4 for Nf = 3, which calls into 
question the notion of an expansion in powers of at low scales. One might also prefer to have 
an expression with free parameters, other than the overall momentum scale A, that could be 
fitted to experimental data. The expressions suggested in Ref. [7] do have free parameters, but 
they also have non-analytic terms, unwanted power corrections and/or singularities at unphysical 
values of k"^. 

2 Model for the running coupling 

If we require 

• No power corrections larger than l/k'^^, 

• Singularities on the negative real k^ axis only, 
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• Some freedom to adjust the form and value at low k"^, 
then the following generalization of Eq. (1.4) suggests itself: 
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1 + c 

\nx ' (l-x)(l + 6) \x + c 
For definiteness we shall consider here 6 = 1/4, c = 4, p = 4, i.e. 



(2.1) 
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(2.2) 



which leads to the form of ag shown in Fig. 1, where the values iV/ = 3, A = 0.25 GeV have 
been chosen, giving as{Q'^) = 0.118 at Q = 91 GeV. For this choice of parameters, as{k^) lies 
very close to the perturbative expression (dot-dashed) above k = 1 GeV, the difference falling 
like Below 1 GeV, however, it remains finite and positive, reaching a maximum value 

of 0.8 at A; = 0.4 GeV. As discussed below, the moments of this function are consistent with 
those deduced from data on power corrections. The multipole on the negative real axis looks 
somewhat unphysical, but should be regarded only as an approximation to a more complicated 
singularity structure. 
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Figure 1: Model for the running coupling (solid curve), compared with one-loop perturbative form 
(dot-dashed) and expansions to first, second and third order in as{Q^) (dashed, Q = 91 GeV). The 
dotted curve shows the 2-loop modification (2.3). 

In order to relate the parameter A in Eq. (1.1) to the scale parameter in some particular 
renormalization scheme, Eq. (2.1) should be generalized to reproduce the 2-loop behaviour of 
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as{k'^) at large k"^. As has been discussed for the simpler form (1.4) [8], this can be achieved 
without spoiling the analytic properties through a replacement of the form 



1 + 



Inx^ ^ 



where 



3(153 - imf) 



(2.3) 



(2.4) 



(33 - 2Nf Y ■ 

As shown by the dotted curve in Fig. 1, this makes only a very small difference at fixed values 
of the parameters, although of course these values would also need some readjustment. For 
simplicity we shall therefore use the 1-loop formula (2.2) in the remainder of this paper. 



3 Power corrections 

As a model for a QCD observable which is of perturbative order a| and is expected to have a 
leading power correction of order 1/Q^, consider the integral [11] 



as{k' 



(3.1) 



Then to lowest order in as{Q'^) we have Fp^q = Fp^q where 



F^)iQ)=\as{Q') 



(3.2) 



Using the explicit expression given by Eqs. (1.1) and (2.2), we can evaluate the integral exactly, 
as shown for p = g = 1 by the solid curve in Fig. 2. The difference between Fp^g and 
arises from the running of as{k'^) in the integrand, which receives both perturbative and non- 
perturbative contributions. It is not possible to disentangle these unambiguously, because the 
two terms in Eq. (2.2) separately give divergent contributions. A perturbative expansion in 
powers of as(Q^) sees only the first term; expanding to 0(a^+^) means using 



asik"^] 
in Eq. (3.1) to obtain 
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which diverges as N oo. This is the infrared renormalon problem [11]-[16]. As one would ex- 
pect from the corresponding curves in Fig. 1, the perturbative estimate (3.4) of Fp^q undershoots 
for low values of N, but becomes arbitrarily large for sufficiently high N (Fig. 2). 
To correct the 0(af+^) perturbative estimate (3.4) we clearly have to write 



Fp,q{Q) = Fi^q\Q) + 6F(^g){Q) 



where the correction term is 



as{k'^ 



a 



(3.5) 



(3.6) 
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Figure 2: Exact vahic of the integral (3.1) for p = q = 1 (solid curve), compared with perturbative 
estimates (3.4) for iV = 0, . . . , 7 (dashed). 



In the approximation scheme proposed in Ref. [2] , one notes that the two terms in the 
of Eq. (3.6) must become similar at values of k in the perturbative region, say k > jXi 
this region cannot contribute to the renormalon divergence. Therefore we neglect the 
above /Xi, the 'infrared matching' scale, and write 



integrand 
, and that 
integrand 



MI 
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where 



N 



(n + g- 1)! 

r being the incomplete gamma function, 

T{n + l,z) -- 
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T{n+l,p\n{Q/ix,)) , 



r e"* dt 
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(3.^ 



(3.9) 



The point of this exercise is that all the necessary non-perturbative information is now contained 
in the Q- and iV-independent parameters Fp^q{iJ,i), which represent weighted averages of aj over 
the infrared region < k < ^i. The resulting contribution to Fp^q{Q) is a power correction 
of order l/Q^- The second term on the right-hand side of Eq. (3.7) subtracts the divergent 
renormalon part of the perturbative prediction. After this subtraction, the total correction 
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q 


p=l 


p = 2 


p = 3 


p = 4 


1 


0.511 


0.450 


0.410 


0.388 


2 


0.283 


0.218 


0.176 


0.155 


3 


0.170 


0.114 


0.081 


0.064 



Table 1: Values of J],,, (2 GeV). 



6Fp^\q) is not simply power-behaved, as may be seen from the difference between the exact 
result and any of the fixed-order curves in Fig. 2. 

For the model of the running coupling shown in 
Fig. 1, we see that Hi = 2 GeV is a reasonable matching 
scale, for which one finds the values of Fp^g{n I J given m 
Table 1. The results Fi,i(2 GeV) ~ ^2,1 (2 GeV) ~ 0.5 
are consistent with the values deduced [1, 2, 3, 5, 14] 
from experimental data on event shapes [17, 18, 19] 
and structure functions [20, 21], respectively.^ 
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Figure 3: Exact value of the integral (3.1) for p = g = 1 (solid curves), compared with perturbative 
(dashed) and infrared matching (dot-dashed) estimates for A'' = 0, 1, 2, 6. 

Using these values in Eq. (3.7), wc obtain the (^Uj = 2 GeV) 'matching estimates' of the 
integral (3.1), shown by the dot-dashed curves in Fig. 3 for iV = 0, 1, 2, 6. The A'^ = 1 case 
corresponds to the most common situation, where we have to match to a NLO perturbative 
calculation. There is still some discrepancy at large Q because the matching in the interval 
IJ,i < k < Q is not perfect. The discrepancy is reduced when the NNLO term is added {N = 2). 
The N = 6 case is shown to illustrate the cancellation of the renormalon and the improvement 
^Note that in Refs. [2, 3] the notation -Fp,i(/ii) = Q;p-i(/ii) was used. 
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in the matching when very many perturbative terms are included. 

Since we know the exact form of the integrand, we could of course improve the matching 
estimates at large Q by choosing a larger value of the matching scale ^i. However in reality 
we usually know only the behaviour as k"^ /Q"^ — > 0, and Hi must be kept small in order for this 
behaviour to predominate in the correction term (3.6). 



4 Renormalons 



A common way to deal with the renormalon divergence of the perturbative estimate (3.4) as 
AT —> DO is to Borel transform and define the sum by the principal value of the Borel integral 
[13, 16, 22, 23, 24, 25]. Recah that if 



then the Borel transform is 



and the sum would be defined as 



F{as) = ^ c„a 

n=0 



n+1 



n=0 



n! 



/•oo 

F{as) = P / duF{u)e-''/'^' 
Jo 



where P indicates the principal value. In the case of Eq. (3.4) we have 

- 1)! P - 2bou 
and hence the 'renormalon estimate' of Fp^q{Q) is 



p 



26oa.(Q2) 



where 



iq-l)\ 

and Ei represents the exponential integral function. 



q-2 
n=0 



/oo p—t 
— dt . 
-z t 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



The uncertainty in the estimate (4.5) is normally taken to be given by the residue of the 
renormalon pole at -u = p/2bo: 



where 



<5Fj-")(Q)=[a,(Q2) 
SHg{z 



P 



2boasiQ^ 



(4.8) 



(4.9) 



(g-1)! • 

For q = 1, the renormalon estimate (4.5) corresponds directly to a principal-value definition 
of the integral (3.1), i.e. to using the perturbative form (1.3) for as{k^) everywhere apart from 
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smoothing the Landau pole over a small region around k = A. This means that as{k'^) < for 
A; < A, leading to a slight underestimation of g when p = q = 1 (Fig. 4). As p increases, 
the negative contribution from /c < A is suppressed and the region of positive dominates, 
leading to a slight overestimate. Nevertheless, for the simple model (2.2) assumed here for as, 
the renormalon estimate is quite satisfactory when q = 1, within the assumed uncertainty (4.8), 
shown by the error bars. 
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Figure 4: Exact values of the integral (3.1) for (solid curves), and renormalon estimates (points). For 
comparison the LO perturbative (dashed) and NLO matched (dot-dashed) estimates are also shown. 



For higher orders, q> 1, the situation is different. If we use the perturbative expression for 
as{k'^), the integrand in (3.1) has a multiple pole and the meaning of the renormalon estimate 
(4.5) is not obvious. A useful way to interpret it (valid also for g = 1) is as follows. The function 
F{Q) = Fp^q{Q) satisfies the differential equation 



p dQ 



asik^) 



which has the general solution 



(4.10) 



(4.11) 
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where Qo is an arbitrary constant of integration. The renormalon estimate (4.5) corresponds to 
a particular choice of Qo > 0) namely the solution of 



z = Zq being the point where Hq{z) = 0, i.e. where 

q-2 

Ei(z) = ^ n! z"'^"^ . 

n=0 



(4.12) 



Explicit values are 



zi = 0.373 , Z2 = 1.347 , z^ = 2.342 . 



0.67 . 



(4.13) 



(4.14) 



Clearly the solution (4.11) corresponds (for Q > 
Qq) to setting asik"^) = throughout < k < Qq. 
Thus the renormalon estimate (4.5) is equivalent to us- 
ing the perturbative expression for as{k'^) everywhere 
above an 'effective cutoff' Qo{p,q) given by Eq. (4.12), 
and setting as{k'^) = below that value. Using the 
1-loop expression for q;s(Qq), we obtain 



q 


p = 1 


p = 2 


p = 3 


p = A 


1 


1.45 


1.20 


1.13 


1.10 


2 


3.85 


1.96 


1.57 


1.40 


3 


10.40 


3.22 


2.18 


1.80 



Table 2: Values of Qo(p,«)/A. 



Qoip,q) = Ae'-'/P ~ /^e^g-o-a-ryp _ 



(4.15) 



Note that the effective cutoff Qo{p,q) grows rapidly with increasing q, especially for p = 1 
(Table 2). 

When q = 1, setting as(k'^) = for k < Qoip, 1) is equivalent to the principal value definition 
of the integral in (3.1), because the principal value of the integral from to Qo{p, 1) is zero. 
For q > 1, the higher effective cutoff leads to underestimation of the integral (Fig. 4). The 
fact that the cutoff on ag is imposed at a scale that depends on p and q makes the renormalon 
approach look unnatural. Given that the analyticity properties are destroyed by the principal 
value prescription, one might as well choose to set as{k'^) = below some fixed value of Qq, to 
be treated free parameter. 

Another way of defining the divergent sum (3.4) is by truncating the series at its smallest 
term [13], i.e. summing up to the integer nearest to 



N{Q) = 



P 



-q + h 



(4.16) 



and assigning an uncertainty equal to the smallest term. This is valid if the series is an asymp- 
totic expansion. The resulting 'optimal truncation' estimate is again a solution of the differential 
equation (4.10), but the effective cutoff Qo{p,q) is now given by iV(Qo) = 0) i-e- 



P 



26o(?-i; 



(4.17) 



which is just Eq. (4.12) with Zq replaced hy q — ^. In view of Eqs. (4.14), one then expects 
results very similar to those of the renormalon approach, which is indeed the case (Fig. 5). 
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Figure 5: Exact values of the integral (3.1) (solid curves), and optimally truncated estimates (points). 
For comparison the LO perturbative (dashed) and NLO matched (dot-dashed) estimates are also shown. 

5 Conclusions 

The simple model of the running coupling introduced in Sect. 2 has allowed us to study various 
aspects of the estimation of QCD power corrections in some detail. For this purpose we have 
used the toy model observable Fp^q{Q) defined in Eq. (3.1). 

The 'infrared matching' procedure introduced in Ref. [2] gives fairly good results when 
matched with the NLO {N = 1) perturbative estimate, improving rapidly as higher orders are 
included (Fig. 3). In principle it eliminates the renormalon problem completely, at the price of 
introducing new non-perturbative parameters, which are the moments of at low scales (Table 
1). The simple model expresses these moments in terms of its own non-perturbative parameters, 
b, c and p in Eq. (2.1). For the present study the latter were fixed at the values in Eq. (2.2), 
which give agreement with a range of data on power corrections. In future these parameters 
could be varied to achieve a best fit. For serious comparisons with data, the 2-loop modification 
(2.3) should be used, and the varying number of active flavours would also need to be taken 
into account. 

A renormalon analysis of the perturbation series suggests another way to estimate power 
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corrections, either through Borel transformation or by optimal truncation of the series. Although 
this approach seems to involve no non-perturbative parameters, it is equivalent to using the 
perturbative expression for the running coupling above a cutoff Qq and setting it to zero below 
Qo- Since the value of the effective cutoff Qo depends on the details of the integrand, viz. the 
powers of momentum and (Table 2) , this procedure cannot correspond to any universal model 
for ag- In particular, the effective cutoff increases rapidly with the power of cXg, and therefore 
renormalon estimates of higher-order effects will fall below those of any universal model for 
sufficiently high orders. 

Whether a universal model of the running coupling, including non-perturbative terms, can 
be extended to higher orders is of course open to question anyway. For the parameter values 
used in Eq. (2.2), the model couphng is nowhere larger than 0.8 (Fig. 1), and so an expansion 
in powers of as{k'^)/iT does not seem unreasonable. The hope is also that large coefficients due 
to infrared renormalons should not appear in this expansion. 
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